Generalized rotating-wave approximation for arbitrarily large coupling 
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A generalized version of the rotating-wave approximation for the single-mode spin-boson Hamil- 
tonian is presented. It is shown that performing a simple change of basis prior to eliminating the 
off-resonant terms results in a significantly more accurate expression for the energy levels of the 
system. The generalized approximation works for all values of the coupling strength and for a wide 
range of detuning values, and may find applications in solid-state experiments. 
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One of the simplest and most ubiquitous models in 
quantum physics is the single-mode spin-boson model, 
consisting of a two-level system coupled to a quantum 
harmonic oscillator. In quantum optics it describes an 
atom coupled to an electromagnetic field mode [l], Q; 
in condensed matter physics it lies at the heart of the 
Holstein model for electrons coupled to phonon modes of 
a crystal lattice Q. More recently, implementations of 
this model have been achieved in superconductor 0, 
6] and semiconductor Q systems. Still other proposals 
have involved mechanical oscillators 8, 9]. Although the 
model itself is quite simple, it displays a rich variety of 
behaviors, encapsulating many of the unique aspects of 
quantum theory. 

The model Hamiltonian may be written as (l7| 
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Despite decades of study, an analytical solution to this 
equation has not yet been found. A number of approx- 
imations have been developed, each tailored to a par- 
ticular range of parameters. In quantum optics, one of 
the most useful approximations is the rotating-wave ap- 
proximation (RWA), which is based on the assumption 
of near-resonance and relatively weak coupling between 
the two systems 0, E| . 

A generalization of the RWA that extends the range 
of validity to arbitrarily large coupling strengths is pre- 
sented in this paper. The only difference from the ordi- 
nary RWA is that a change of basis is performed prior 
to carrying out the approximation. For the case of exact 
resonance (O = luq), the energy levels given here were 
first found by Amniat-Talab et al. However, their 

derivation involved a complicated method of quantum 
averaging and resonant transformations. The derivation 
presented here is not restricted to exact resonance and 
the resulting approximation works remarkably well for 
large detuning. Moreover, in this form the simplicity of 
the approximation and its close connection to the stan- 
dard RWA are emphasized. 

To begin with, a brief review of the standard RWA is 
given in order to establish the arguments used in deriv- 
ing the generalized approximation. The first step is to 



rewrite Eq. JT]) in the form 

H = LjQaJa + \£lcr x + A(<7_<v + <r + a + & + a) + <r_a) 

(2) 

where a± — \{<t z T i&y) are the raising and lowering 
operators in the basis of a x . Alternatively, the Hamilto- 
nian may be written in matrix form in the basis |±x, AT) 
(where N — 0, 1, 2, ... ), which is the eigenbasis of the 
noninteracting Hamiltonian Hq = ojqo) a 
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where E 



(0) 

±,N 



Ncuo ± ifi and the order of the columns 
and rows is \—x, 0) , \+x, 0) , |— x, 1) , \+x, 1) , — 

Consider the case of near-resonance (u>o ~ ^) and weak 
coupling (A -C u}q,£1). The interaction term <r_at + <r + a 
couples the states \+x,N) and |— x, N + 1), which have 
nearly equal energies in the absence of the interaction. 
On the other hand, the term <j + ci i + er_a couples the off- 
resonant states \—x,N) and \+x,N + 1). In this sense 
the first term is "energy conserving," while the second 
is not. The rotating-wave approximation eliminates the 
non-energy-conserving terms. In matrix form this cor- 
responds to removing the remote matrix elements. The 
Hamiltonian then becomes block diagonal and may be 
readily diagonalized. 

Alternatively, the RWA Hamiltonian may be derived 
by moving to the interaction picture with respect to Hq. 
The Hamiltonian becomes 



H[{t) = exp(ii7 £)A<7 z (a^ + a) exp(~iH t) 



(4) 
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In the case of near resonance, ujq w n and the first two 
terms vary slowly in time. The last two terms, however, 
vary rapidly and therefore average to zero over timescales 
on the order of 1/u)q. Thus the last two terms may be 
neglected and, transforming back to the Schrodinger pic- 



ture, the coupling reduces to A(<7_ 



it 



■ a+a). These are 



exactly the "energy conserving" terms discussed in the 
matrix derivation. Indeed, the argument about energy 
conservation and the argument about timescales are com- 
pletely equivalent. 

The next task is to establish a new set of basis states. 
These are the states obtained from an adiabatic approx- 
imation in the limit Q <C (<x>o,A). Several derivations 
of this approximation have been presented 




However, the matrix-based derivation given in Ref. \U 
is the most useful for the purposes of this paper and is 
briefly summarized here. 

The adiabatic approximation is most readily carried 
out in the basis obtained by setting £1 = in Eq. ([1]): 



\±z,N ± ) = \±z) 
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(5a) 
(5b) 



The qubit states |±z) are eigenstates of a z and the oscil- 
lator states \N±) are position-displaced Fock states. Note 
that \+z, N+) and |— z, 7V_) are degenerate in energy. 

The spin term \^1<J X couples the basis states given in 
Eq. (f5a[) . Within the adiabatic approximation, only the 
coupling between states with the same value of N is con- 
sidered. In matrix form, this corresponds to reducing 
the matrix to a block diagonal form, where the blocks 
are given by 
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The expression (N- \N + ) is simply the overlap of the two 
position-displaced Fock states, given by (M < N) 



The 2x2 matrix of Eq. 
ergies 



(7) 

has the eigenstates and en- 



\*±,n) = jsil+z, N + ) ± |-2, N_)), (8a) 
E±, N =±fa{N-\N+)+E N . (8b) 

An analysis of the adiabatic approximation and its con- 
sequences may be found in Ref. [14j. 

The derivation of the generalized rotating-wave ap- 
proximation (GRWA) is now quite straightforward. The 
Hamiltonian is rewritten in the basis of the adiabatic 
eigenstates |^±,jv)- Then the argument about energy 
conservation that led to the RWA is applied in the new 
basis and the approximate energy levels are calculated. 



When written in the basis of the states 
|*-,o) , |*+,o) , , , ■ ■ ■ , Eq. © becomes 
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where n 



M.N 



(9) 

n(M_\N + ). The form of this matrix 
closely resembles that of Eq. ^ with additional remote 
matrix elements. As before, the approximation consists 
of neglecting the remote matrix elements, reducing the 
matrix to a 2 x 2 block diagonal form. 

Although it is not immediately evident from the ma- 
trix form, the terms retained in this approximation corre- 
spond to energy-conserving one-particle transitions, just 
as in the ordinary RWA. This is most easily illustrated 
in the interaction picture. First the change of basis from 
\^fx,N) to I'l'^jv) is carried out by a unitary transfor- 
mation with the operator 



'St 



(10) 



which is a spin-dependent position displacement opera- 
tor. Applying this transformation to Eq. |T]) results in 
the transformed Hamiltonian 

H = b ] HD 

= w a t a+|^xexp[-^(T 2 (ffi , -a)] (11) 
= ujQa)a + jCla x + Hi. x + Hi, y , (12) 

where the functions H\ x and H\ y are defined as 



H hx = \m x 2(A)* (a t - af + |(^) 4 (at - af + 



2(A )( at_ 6 ) + |(JL)3(at-a) 3 + 



(13) 



(14) 



The next step is to move to the interaction picture with 
respect to Hq = ujqwcl + hQa x , Let us examine H\ }X 
first. Since a x commutes with the rotation operator 
U = cxp(iHot), the rotation affects only the oscillator 
operators. Take the first term of H\ x as an example. In 
the interaction picture the operators become 



d x {— 2a)a — 1 



a t2 e 2^o* 



a 2 e 



2iujQt 



(15) 



The time-independent terms contain powers of the num- 
ber operator a) a and correspond to transitions that result 
in zero net excitation of the oscillator. They are diagonal 
in the basis |$±,j\r) and modify the spin frequency n, re- 
sulting in the term ±in (N_\N + ) that appears in E± t N- 
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The higher-order terms, which have a rapid time depen- 
dence, produce remote matrix elements such as the two- 
excitation term (^_^ \H\^f_ t2 ) and are neglected within 
the GRWA. 

Next consider H\ y . Taking ia y = — &+), the first 
term is proportional to 



(16) 



The second term of Hi, y is slightly more complicated. 
When (a' — a) 3 is expanded and put into normal order, 
this term is given in the interaction picture by 

[ - 3a-tf(tfa + l) e ^°-n)i + ^ + {a)a + l)ae-^°- n)t 
- M + tf(tfa + l)e i(wo+n)t + 3<r_(a t a + l)ae- i( "*> +n)t 

+ ^_ a t3 e <(3<-o-n)t _ - +a 3 e -i(3w -Q) 4 

(17) 

The first two terms create energy-conserving transitions 
involving a single excitation. They produce the matrix 
elements + , N \H\ty -, N+ i) and (*_ iiv +i|i?|*+,Ar) that 
appear immediately off the diagonal in Eq. . The next 
two terms correspond to energy non-conserving single- 
excitation transitions and produce the remote matrix el- 
ements w|if|* +! ./v + i) and (^+M+i\H\^- tN ). The 



last four terms involve a net change of three excitations 
and produce remote matrix elements. Only the first 
two terms of Eq. (jTTJ) have slow time dependence when 
CI « ujq. 

Finally, the GRWA is carried out by keeping only 
the "energy-conserving" one-excitation terms. The other 
one-excitation terms as well as terms involving higher 
numbers of quanta are discarded. When all powers of 
X/luq are taken into account, H\, y reduces to a coupling 
term of the form 



n(£)[<7_at/(ata) + * + /*(fit a )a] 



(18) 



where the function f(a'a) is too complicated to display 
here. Equation l|18j) is a generalization of the energy- 
conserving term \(&_aJ + a+a) in the usual RWA Hamil- 
tonian 

Returning to the matrix picture, the GRWA ground 
state is uncoupled from all the other states (just as in 
the RWA), so the ground state energy is given by -E-,0- 
The remainder of the matrix takes the familiar 2x2 
block-diagonal structure with blocks of the form 



E+,N-1 2^N-l,N 



±ci' 



N 



(19) 



Solving for the eigenvalues of the blocks yields the GRWA 
energies: 



£GRWA = {N + l )wo _ A_ + U _ 2A2/W? [Ln[A}?/uj1) _ iA , +l( 4 A 2 /w 2 )] 

± Qlo, -lne- 2 ^ [L N (4X 2 /u, 2 ) + L N+1 (4\*/J&)] }' + ^ [JtJ r (4A a /«g)] a ) ^ 

(20) 

I 



The energy levels from the RWA, the adiabatic approxi- 
mation, and the GRWA are plotted in Fig. Q] For com- 
parison purposes, the energy levels obtained from a nu- 
merical solution of Eq. (JTJ) are also shown. The RWA 
reproduces the correct limiting behavior as X/loq — > 0, 
but breaks down near the point where the paired levels 
first cross. On the other hand, the adiabatic approxima- 
tion diverges from the numerical solution at small values 
of X/loq, but captures the behavior beyond the first cross- 
ing point very well. The GRWA combines the behavior of 
the adiabatic approximation at large values of A with the 
accuracy of the RWA at small values, providing an excel- 
lent approximation to the actual energies of the system 
over the full range of coupling strengths shown. 

Remarkably, the GRWA works reasonably well even for 



large detunings with ojq < CI. As an example, the case 
u)q — 0.75CI is illustrated in Fig. [2] The maximum error 
in the energy is less than 0.2lo for the ground state and 
decreases for higher energy levels. The qualitative agree- 
ment between the GRWA and the exact solution remains 
fairly good even down to loq = 0.5f2. Considering that the 
RWA requires small detuning and the adiabatic approxi- 
mation is derived under the assumption that luq 3> Cl, the 
GRWA is surprisingly robust in this parameter regime. 

Why does the GRWA work so well? It seems coun- 
terintuitive that a simple change of basis for the RWA 
should result in such an improvement. One explanation 
comes from examining a fundamental similarity between 
the RWA and the adiabatic approximation: both involve 
calculating the energy splitting due to an interaction be- 
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FIG. 1: Comparison of the RWA (dot-dashed), adiabatic ap- 
proximation (dotted), and GRWA (dashed) with numerically- 
determined energy levels (solid) in the resonance case, u>o = 

n. 



are capable of much higher coupling strengths, even at 
large detunings, and are already nearing the limits of va- 
lidity of the RWA dill. 

The generalized approxima- 
tion presented in this paper provides an excellent treat- 
ment for the very strong coupling limit that these exper- 
iments are expected to achieve, while maintaining strong 
links to the familiar language and techniques of quantum 
optics. Thus the GRWA may prove useful as experiments 
continue to expand the accessible parameter regime in 
this important and still fascinating model. 
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M. S. Kim, and M. Paternostro for their helpful com- 
ments and encouragement. Support from the National 
Physical Sciences Consortium and the EPSRC is grate- 
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FIG. 2: GRWA energy levels (dashed lines) compared 
with numerically-determined energies (solid lines) in the off- 
resonance case, loq = 0.75J1. 



tween two otherwise degenerate basis states. In the RWA, 
the degeneracy or resonance occurs at the single point 
(f2 = u>q, A = 0). The adiabatic approximation, on the 
other hand, treats the resonance at Q = 0, which occurs 
for all values of A. This interpretation accounts for the 
fact that the RWA only works for small A as well as for 
the accuracy of the adiabatic approximation at all values 
of A when fl <C luq. The adiabatic approximation breaks 
down at small A when Q = loq precisely because it does 
not account for the zero-coupling resonance [h|. How- 
ever, the GRWA takes into account both the resonance 
at £1 = and the point-like resonance at (il = a->o, A = 0), 
which yields a very accurate energy spectrum. 

One reason the standard RWA has remained so preva- 
lent in quantum optics is that the conditions of near- 
resonance and weak coupling are naturally satisfied in 
atomic cavity quantum electrodynamics (CQED) exper- 
iments 
systems 
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